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We present a comparative theoretical study of the effects of standard Anderson and magnetic 
disorders on the topological phases of two-dimensional Rashba spin-orbit coupled superconductors, 
with the initial state to be either topologically trivial or nontrivial. Using the self-consistent Born 
approximation approach, we show that the presence of Anderson disorders will drive a topological 
superconductor into a topologically trivial superconductor in the weak coupling limit. Even more 
strikingly, a topologically trivial superconductor can be driven into a topological superconductor 
upon diluted doping of independent magnetic disorders, which gradually narrow, close, and reopen 
the quasi-particle gap in a nontrivial manner. These topological phase transitions are distinctly 
characterized by the changes in the corresponding topological invariants. The central findings made 
here are also confirmed using a complementary numerical approach by solving the Bogoliubov-de 
Gennes equations self-consistently within a tight-binding model. The present study offers appealing 
new schemes for potential experimental realization of topological superconductors. 
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Introduction. —Topological superconductors (TSCs) 
[1] have been intensively explored recently as candidate 
systems for realization of Majorana fermions [2]. Due to 
their exotic non-Abelian braiding statistics [3l|4], Majo¬ 
rana fermions are in turn expected to play an important 
role in future fault-tolerant topological quantum compu¬ 
tation HIE]. Many different schemes have been proposed 
to realize TSCs, including centrosymmetric odd-parity 
pairing copper oxide superconductors miHi, surface states 
of topological insulators laiini or two-dimensional (2D) 
Rashba spin-orbit coupled semiconductors m proximity 
coupled with an 5-wave superconductor, and ID ferro¬ 
magnetic Shiba chain on top of a conventional supercon¬ 
ductor with strong spin-orbit coupling [12]. These dif¬ 
ferent innovative proposals continue to stimulate active 
research efforts on definitive experimental realization of 
TSCs and observation of Majorana fermions. 

In essentially all the compelling experimental demon¬ 
strations of TSCs and/or Majorana fermions reported so 
far [T 2 UT 6 ] , the presence of certain types of disorders had 
to be unavoidable. Generally speaking, disorders are un¬ 
desirable physical entities, because they may destroy the 
salient properties of a clean system. On the other hand, 
there have also been numerous examples showing that 
properly introduced disorders can lead to the emergence 
of various intriguing new physical phenomena in other¬ 
wise clean but ordinary host systems. One such example 
is the recent discovery of the topological Anderson in¬ 
sulating (TAI) state in HgTe quantum wells [IllllH]. It 
was shown that the presence of Anderson disorders can 


drive a topologically trivial insulating state into a topo¬ 
logically nontrivial state (or the TAI state), with a quan¬ 
tized conductance of Go = jh [T7II2Q] . In view of the 
widespread presence of disorders in realistic systems and 
the recent example on TAI, it is natural to investigate the 
feasibilities of disorder-induced topological phase transi¬ 
tions in superconducting systems, especially the potential 
disorder-assisted realization of TSCs. 

In this Letter, we carry out a comparative study of 
the effects of standard Anderson and magnetic disor¬ 
ders on the topological phases of 2D Rashba spin-orbit 
coupled superconductors using complementary theoreti¬ 
cal approaches, with the initial state to be either topo¬ 
logically trivial or nontrivial. First, we use the self- 
consistent Born approximation (SCBA) to investigate 
the disorder renormalized density of states (DOS), and 
show that the presence of usual Anderson disorders will 
drive a topological superconductor into a topologically 
trivial superconductor in the weak coupling limit. More 
strikingly, a topologically trivial superconductor can be 
driven into a topological superconductor upon doping 
of diluted magnetic disorders, which drive an intricate 
closing-and-reopening process of the quasi-particle gap 
in a nontrivial manner. These topological phase tran¬ 
sitions are quantitatively characterized by their corre¬ 
sponding topological invariants, and the central findings 
are also confirmed by solving the Bogoliubov-de Gennes 
equations self-consistently within a tight-binding model. 
The present study offers new insights towards potential 
experimental realization of TSCs. 
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Theoretical model —To start, we consider a 2D elec¬ 
tron gas with a Rashba-type SOC and a Zeeman field 
h. Discussions on potential experimental realizations of 
such systems will be deferred to near the end of this pa¬ 
per. By introducing a proper attractive interaction be¬ 
tween the electrons, the system is in a superconducting 
state, and can further be classified as a topological su¬ 
perconductor when the Zeeman field exceeds a critical 
value he mumna. The total Hamiltonian H of such 
a 2D system contains two parts. The first part is the ki¬ 
netic energy, given as Hq = f dr'ip^r)'Ho(r)'i/j(r), where 
V’^(r) = (V’|(r),'*/)j(r)) and ?^o(r) = j{2me) - 

i\{dy(7x — dxCFy) — h(7z. Here, rUg is the effective mass 
of the electrons, fi is the chemical potential, A is the 
strength of the Rashba SOC, and with i = x^y^z are 
the Pauli spin matrices. The second part is the on-site 
attractive interaction between the electrons, described as 
Ha = U J dr'0|(r)'0|(r)'0^(r)'0^(r), where U{< 0) is the 
attraction strength. In the momentum space, by treat¬ 
ing the two-body interaction Ha within the mean-field 
approximation, the total Hamiltonian H can be written 
as = 1/2 Y.k with 

%{%) = {€k - + \{kyCJa; “ k^CTy) - ha^ + AcT^Ta;, (1) 

where = (c^^,are the field opera- 

tors in the Nambu spinor basis, Ck = j2mn with 

i = x,y, z are the Pauli matrices acting on the particle- 
hole degrees of freedom, and A = U ^j^{c-kiCk^) is the 
mean-field superconducting order parameter. The sys¬ 
tem described by Eq. (1) is a topologically non-trivial 
superconductor when h > he = [23]. Here 

we differentiate the present model systems, where time- 
reversal symmetry is explicitly broken by the Zeeman 
term, with those considered in some earlier studies in 
which time-reversal symmetry is protected [24l . 

To describe the interactions between the electrons and 
disorders introduced into the system, we consider a lo¬ 
cal scattering Hamiltonian Him = — f dr'lp^r)Vimd(r — 
ro)'0(r), whose simple form is able to capture the essen¬ 
tial physics to be exploited in the present study. In the 
case of the usual Anderson disorders. Vim is the on-site 
disorder potential V distributed uniformly in the interval 
(—Vb, Vb). For magnetic disorders, we choose an isotropic 
form given by Vim = J{S • a), where J denotes the ex¬ 
change coupling strength between the electrons and mag¬ 
netic disorders, a = (da,, (Ty^dz) for the electron spin, and 
S is the moment of a magnetic dopant. In Nambu nota¬ 
tions, the disorder scattering Hamiltonian can be rewrit¬ 
ten as 

-ffim = -^ E (2) 

^ kk' 

where Vim = for the Anderson disorders and Vim = 
JS -a for the magnetic disorders. Here the effective elec¬ 
tron spin operator a is defined as a = l/2[(l-\-Tz)d-\-{l — 


Hamiltonian (2) is a general form of the inter¬ 
action between the disorders and superconducting quasi¬ 
particles, and Eqs. (1) and (2) provide the framework for 
further studying the disorder effects on the topological 
phases of the superconducting systems. 

The Green’s function formalism is well suited to study 
the effects induced by multiple disorders in a supercon¬ 
ductor. After averaging over the randomly distributed 
disorders, the Matsubara Green’s function of the system 
described by Eqs. (1) and (2) is given as 

G-\k,UJn)=i^n-'Hik)-t, (3) 

where uJn = (2n+l)7r/cBT, ks is the Boltzmann constant, 
T represents the temperature, and E is the self-energy. 
In the presence of disorders, the superconducting order 
parameter is determined by the self-consistent equation 

A = -if/fcBT^^lV[G(fc,w„KrJ. (4) 

k 

In the SGBA approach, E is given by 

E(A^,CJ^) "klim^^VimGik^ UJu^Vim-) (5) 

k 

where riim is the concentration of the disorders. Taking 
into account of the symmetry restriction and the matrix 
structure of E, the self-energy effects can be manifested 
as disorder-induced renormalizations of cjn, hi and A. 
As shown in Eq. (5), E is independent of the momentum, 
leading to a renormalized form of the Green’s function 
G{k,Un) = [ioJn - (e/c - h)rz - \{ky(Jx - kxCJy) + hcTz - 
obtained with the replacements of ^ 
y ^ h ^ h, and A ^ A. The solution of Eq. (5) 
gives rise to a set of self-consistent equations 

UJn = UJn aiFi (uJni h, A^ 
h = h - (^2^2 (^UJni hi A^ h, (6) 

A = A + q^s-Fs (^uJni h^ A^ A, 

where 0 ^ 1 , 2,3 = -^nimVQ for the Anderson disorders, and 
Qfi = -0^3 = \nimJ‘^S‘^, a2 = \nimJ‘^{Sl F si - S^) 
for the magnetic disorders. Here we have derived the de¬ 
tailed analytical formulas for Fi^ 2 , 3 (^n, A) from Eq. (5) 
by considering the matrix structure of the Green’s func¬ 
tion (given in the Supplemental Material [26]). We also 
ignore the renormalization of jd as it simply changes the 
chemical potential, which can be preset suitably, unlike 
the TAI cases, where the renormalization of jd has signifi¬ 
cant effects [T8ll2Q] . The disorder-averaged Green’s func¬ 
tion G{k,uJn) can then be obtained by solving Eqs. (6) 
and (4) self-consistently. 

Topological phase transitions. —To study the disor¬ 
der effects, we investigate the DOS defined as p{uj) = 
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FIG. 1: (Color online) Topological phase transitions induced 
by the Anderson disorders, with the initial TSC state obtained 
by choosing h = 3.5 meV. (a) The DOS as a function of the 
disorder strength Vq. The dashed (green) curves highlight 
the edges of the quasi-particle gap. (b) Topological invariant 
as a function of Vo at different A, with the solid (red) curve 
corresponding to the phase transition shown in (a). Other 
parameters include rrie = 0.01 meV-A“^, and U = 160 meV, 
leading to A 3 meV. 


— (l/27r)Im{Tr[Gr(co’)croro]}, where Gr(co’) is the retarded 
Green’s function obtained by performing analytical con¬ 
tinuation on G{uJn)- According to Anderson’s theorem 
m, a conventional 5-wave superconductor can be barely 
influenced by potential disorders (later known as Ander¬ 
son disorders), but such an impurity has been shown 
more recently to be able to induce mid-gap bound states 
in a TSC [28l[29]. It is therefore worthwhile to study the 
effects of multiple Anderson disorders in TSCs relative 
to that in normal superconductors. Furthermore, a mag¬ 
netic impurity can lead to the appearance of the Yu-Shiba 
states |30l [31] inside the superconducting gap even if the 
superconducting state is topologically trivial, making it 
more appealing to study the effects of multiple magnetic 
impurities in different classes of superconductors. 

Here we first investigate the effects of multiple Ander¬ 
son disorders on the topological phases in TSCs. Based 
on the study of the DOS, we recover a topological phase 
transition from the TSC state to a trivial superconduct¬ 
ing state as illustrated in Fig. 1 (a). Qualitatively, an 
increase in the strength Vb of the on-site Anderson dis¬ 
orders will gradually close the quasi-particle gap of the 
TSC. When Vb crosses a critical value Vb, the quasi¬ 
particle gap reopens, and the system adopts a trivial 
superconducting state. These findings suggest that, as 
far as Anderson disorders are concerned, a cleaner su¬ 
perconducting system is better for potential realization 
of TSC. Further increment of Vb to the strong disorder 
regime is believed to transform the system to the Ander¬ 
son localized state [32] . 

We next study the effects of the magnetic impurities. 
Unlike the case of the Anderson disorders, the presence of 
the magnetic disorders in a TSC cannot drive any topo¬ 


FIG. 2: (Color online) Topological phase transitions induced 
by the magnetic disorders, with the initial trivial SC state ob¬ 
tained by choosing h = 2.5 meV. (a) The DOS as a function of 
Him with J = 100 meV. The dashed (green) curves highlight 
the edges of the quasi-particle gap. (b) Topological invariant 
as a function of riim at different A, with the solid (red) curve 
corresponding to the transition shown in (a). Other parame¬ 
ters include S = v^, rUe = 0.01 meV-A“^, and U = 160 meV, 
leading to A 3 meV. 

logical phase transitions, indicating that the TSC state is 
relatively robust against magnetic disorders. Remarkably 
and most strikingly, we find that an initial topologically 
trivial superconductor can be converted into a TSC via 
diluted doping of randomly distributed magnetic impu¬ 
rities. As shown in Fig. 2 (a), the system undergoes a 
gap closing and reopening process upon increasing the 
concentration riim of the magnetic disorders. A critical 
concentration Uc for such a quantum phase transition 
is clearly identifiable. Another feature of Fig. 2 (a) is 
that the quasi-particle gap closes again upon further in¬ 
creasing riim- In this gapless region containing sufficient 
magnetic impurities (with Uim > 1 . 2 % in the present 
study), the system is a gapless superconductor with non¬ 
vanishing superconducting order parameter [33]. How¬ 
ever, below this concentration regime, we do have a siz¬ 
able range of during which TSCs can be observed. 

To measure the occurrence of these topological phase 
transitions induced by disorders, we further calculate the 
topological invariant for various parameter regions. Due 
to the presence of the disorders, their interaction with 
the host system makes it impossible to use the stan¬ 
dard Berry phase approach based on a pure band picture 
in calculating the Chern number [34]. Here we adopt 
an alternative and more general formula for evaluating 
the topological invariant, which relies on the full Green’s 
function of the interacting system as [354(37] : 

A/b = — J Tr[eQ,/3^G9Q,G ^Gdj^G ^Gd^G ^], (7) 

where 60 ,^ 3 ^ is the Levi-Civita symbol, p = (cj^, /c), dcx = 
dp ^, and the summations over [3^^ are implied. The 
full Green’s function of such a disordered system can be 
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obtained self-consistently, thus the topological invariant 
is a straightforward calculation using Eq. (7). As shown 
in Fig. 1(b), the topological invariant jumps from A/2 = 1 
to 0 at the critical value of Vb, which indicates a topolog¬ 
ical phase transition from a chiral TSC to a trivial SC. 
Figure 2(b) further confirms that a topologically trivial 
SC (A/2 = 0) can be converted to a chiral TSC (A/2 = 1) 
via increasing magnetic doping. In addition, we also find 
that the enhancement of the SOC strength will promote 
these topological phase transitions in both cases of the 
Anderson and magnetic disorders. Collectively, these re¬ 
sults offer strong evidence for the existence of rich topo¬ 
logical states induced by proper choices of disorders. 

Numerical BdG solutions .—In this section, we numeri¬ 
cally investigate the disorder-induced effects within a cor¬ 
responding tight-binding model of the system described 
by Eq. (I). By projecting Hamiltonian (I) on a 2D square 
lattice, we reach 


H = -tJ2 (clcj. + H.c.) + A + H.c.) 

(b')cr {ij) 

h ^ ^ U ^ ^ ^ ^ '^icr 

i(j<j' i ia 

T ^ crcr'C'ifjC'icr' •! 

( 8 ) 

where t is the nearest-neighbor hopping term, ct^ and 
Qcr are the creation and destruction operators for an 
electron with spin a on site of the square lattice, 
h-icr = ^ia^i(T is the particle number operator, and is 
the angle between (r^ — r^) and the x axis. Here, we solve 
the disordered system characterized by Eq. (8) within 
the standard mean-field BdG approach, 7^BdG^n(i*i) = 
where 


^BdG = 


0 




7^ 0 

niiYi) -A(ri) 
-A*(r,) -Ulin) 
0 


\ 

0 


(9) 


is the BdG Hamiltonian, is the eigenenergy of the 
corresponding quasi-particle wave function T^(r^) = 
(r^), (r^), z^I(r^), z/I(r^)]^ defined in the Nambu 

spinor space. Here 'H^{Yi)u^{Yi) = —t^^u^{Yi -|- (5) + 
7^(rjX(rj) = + 5)- 

S = ±x,±y, A(rj) is the local superconduc- 
tor order parameter, and /i^o- = /a-\-U{ni,-a) is the mod¬ 
ified chemical potential due to Hartree shift. We note 
that = J{Sx + iSy) for the magnetic disorders and 
vanishes for the Anderson disorders. Starting with some 
initial guess values of A(r^), we first numerically solve 
the BdG Hamiltonian on a A x A" square lattice with pe¬ 
riodic boundary conditions. Next, we calculate the local 
pairing amplitudes and particle density given by 


M^i) = U riia 

n 


( 10 ) 



FIG. 3: (Color online) The upper and lower panels show the 
numerical results of a 60 x 60 square lattice in the presence 
of the on-site Anderson and magnetic disorders, respectively, 
(a) Quasi-particle gap as a function of the disorder strength 
Vq. (b) The spatial variations of A(rA for the Anderson dis¬ 
orders with Vb = 10 meV. (c) Quasi-particle gap as a function 
of the on-site magnetic coupling strength J. (d) The spatial 
variations of A(ri) for randomly distributed magnetic disor¬ 
ders with J = 4.8 meV. The error bars in (a) and (c) show 
the standard deviations of the mean for 10 samples, and the 
inserts give the averaged local superconducting order param¬ 
eters. Other parameters include S = a/6, U = 160 meV, and 
t = 50 meV, corresponding to me — O.OI meV-A“^. 


and iterate this process until self-consistent values of 
and A(r^) at each site are achieved. 

Numerical solutions of the BdG equation can yield the 
whole quasi-particle spectrum in the presence of the dis¬ 
orders. As shown in Figs. 3 (a) and (c), there exists two 
regimes of the disorder strength showing nonzero quasi¬ 
particle gaps. Similar to the SGBA results, we also ob¬ 
serve gap closing and reopening processes by increasing 
the strength of the disorders. These behaviors are quali¬ 
tatively the same as that shown in Fig. 1(a) and Fig. 2(a). 
From the quantitative perspective, the critical value of 
Vo for the Anderson disorders is about 26 meV, which 
is comparable to Vb = 22 meV shown in Fig. 1(a). For 
computational convenience, the magnetic disorders are 
treated as on-site spins with randomly oriented directions 
in the numerical calculations. As shown in Fig. 2(a), 
the critical concentration characterizing the topological 
phase transition is Uc = 0.4%, translating into J = 3.2 
meV for on-site disorders, and the value is comparable 
to the critical coupling strength = 3.6 meV illustrated 
in Fig. 2(c). Based on these qualitative and quantita¬ 
tive comparisons between the results obtained analyti¬ 
cally earlier and numerically now, we can conclude af- 
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firmatively that rich topological states can indeed be in¬ 
duced by properly introducing disorders. In particular, 
a topologically trivial SC can be readily converted into a 
TSC upon diluted doping of independent magnetic impu¬ 
rities. In addition, the spatial variations of A(r^) shown 
in Figs. 3(b) and (d) demonstrate that the systems evolve 
into phase-separated regions in the presence of disorders, 
and the sizes of these regions are comparable to the su¬ 
perconducting coherence length ^ = up/A ^ 10a, where 
uf is the Fermi velocity, and a is the lattice constant. 

Before closing, we briefly discuss candidate systems 
for potential experimental realizations of such disorder- 
induced topological phase transitions. We could consider 
a superconductor thin film with strong SOC, such as Pb 
[121 EHl [39] or PbBi alloyed films grown on semiconduct¬ 
ing substrates. Due to the lacking of the inversion sym¬ 
metry, the Rashba-type SOC will be present on the sur¬ 
face of such superconductors. Therefore, by directly dop¬ 
ing magnetic elements into or on the surface of such 2D 
superconductors, their topologically trivial nature may 
be converted into TSCs. Our results may also be ob¬ 
served in the recently realized SOC coupled 2D ultracold 
atomic Fermi gases of [40] and ^Li atoms m, which 
are shown to be topological superfluids. In addition, the 
present findings can be extended to ID and 3D cases. 

In summary, we have studied, using complementary 
analytical and numerical approaches, the effects of the 
Anderson and magnetic disorders on the topological 
phases of 2D superconductors with Rashba SOC. We 
have found that the presence of the Anderson disorders 
will drive a TSC into a topologically trivial SC in the 
weak coupling limit. More strikingly, a topologically triv¬ 
ial SC can be converted into a TSC upon diluted doping 
of independent magnetic impurities, which is character¬ 
ized by an intricate nontrivial gap closing and reopening 
process. The central findings can offer new insights to¬ 
wards potential experimental realization of TSCs. 
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